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Abstract 

Within a QCD-based eikonal model with a dynamical infrared gluon mass scale we discuss how 
the small x behavior of the gluon distribution function at moderate is directly related to the 
rise of total hadronic cross sections. In this model the rise of total cross sections is driven by 
gluon-gluon semihard scattering processes, where the behavior of the small x gluon distribution 
function exhibits the power law xg{x, Q'^) = h{Q'^)x~''. Assuming that the scale is proportional 
to the dynamical gluon mass one, we show that the values of /i(Q^) obtained in this model are 
compatible with an earlier result based on a specific nonperturbative Pomeron model. We discuss 
the implications of this picture for the behavior of input valence-like gluon distributions at low 
resolution scales. 
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I. INTRODUCTION 



The increase of hadronic total cross sections was theoretically predicted many years ago 
|l| and this prediction has been accurately verified by experiment I2]. Nowadays, one of 
the main theoretical tools to explain this behavior is the diffractive formalism underlying 
the so called QCD-inspired models j^, 0, [s, [g]. In this formalism the scattering amplitude 
in the cm. system is associated with semihard processes, i.e. processes with very small x 
whose amplitudes can be calculated in perturbative QCD: at high energies the typical trans- 
verse momentum of semihard processes is relatively large and for wee partons in hadronic 
scattering processes it leads to small values for the strong coupling constant ag- At high en- 
ergies semihard processes compete with "soft" ones which have traditionally been regarded 
as being responsible for the bulk of the hadronic cross section. Hence the study of hadronic 
scattering processes is no further away from the QCD theory than, for example, the study of 
hard processes at small transverse distances; as the energy increases, semihard processes are 
expected to give an increasing and significant part of the total hadronic cross section 0, Q] . 
This appealing idea is a feasible one since it is well known that experimental observations 
made it clear that at high energies the soft and the semihard components of hadron-hadron 
scattering are closely related, just as, for example, the observation of correlations between 
the average transverse momenta of hadrons produced and the multiplicity density in rapid- 
ity, SI and the copious production of gluon jets with moderate pt (minijet phenomena) in 
hadronic collisions [lol |. 

QCD-inspired models incorporate soft and semihard processes in the treatment of high 
energy hadron-hadron interactions using a diffraction-scattering formulation compatible with 
analyticity and unitarity constraints. At large energies and fixed transverse momentum there 
is a rapid growing number of semihard gluons in the hadron, which is the main reason of 
the rising total cross section. Therefore the rise of total hadronic cross section is intimately 
connected to the small x behavior of the gluon distribution function and, within the range 
of applicability of these models, it can be used to investigate the gluon distribution exactly 
in the x <C 1 semihard region. 

In QCD-inspired models a phenomenological gluon distribution with the behavior g{x) ~ 
hx~^ is used in the actual calculations, where J > 1. It could be asked why this approach 
has not been used before to determine the small x behavior of g{x). One possible reason 
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lies in the fact that previous models were quite dependent on an ad hoc infrared cutoff mass 
scale and on the infrared value of the strong coupling constant (ao) Therefore, 
if g{x) at small x is given by the above expression, the determination of h would not be 
reliable since it appears multiplied by ao, which is a fitting parameter in the model. 

The freedom to choose the values of the mass scale and disappears in one iniprovement 
of the model, where the arbitrary mass scale is changed by a dynamical one (6|. In this 
model (henceforth referred to as DGM model) the onset of the dominance of semihard 
gluons in the interaction of high-energy hadrons is managed by the dynamical gluon mass 
rrig [Ul (intrinsically related to an infrared finite gluon propagator, |l2l |) whose existence is 
strongly supported by recent QCD lattice simulations [l3| as well as by phenomenological 
results for recent reviews on the phenomenology of massive gluons 



see Ref. [18| and references therein. In the DGM model the infrared coupling constant 
is also a function of the dynamical gluon mass Therefore we have a precise physical 
meaning for this infrared coupling as well as for the quoted dynamical mass scale, which is 
a natural regulator of the infrared divergences associated with the semihard gluon-gluon 
subprocess cross sections. These properties render the DGM model instrumental for a 
reliable determination of the QCD-predicted small x distribution gi^x^Q"^) = h{Q'^)x~'' at 
moderate Q^, namely the determination of the h{Q'^) factor, where the Q scale is on the order 
of the dynamical gluon mass m^. This dynamical scale represents the onset of semihard 
contributions to hadronic diffractive scattering and hence provides the necessary scale in 
order to apply the diffractive formalism. 

The gluon distribution function is usually determined in the large region and its small 
X limit in this region has been extensively discussed in recent times. Therefore the results 
that we will discuss in this work can be considered phenomenologically important if we want 
to understand the match between these different kinematical regions. Moreover, the gluon 
distribution at a low value of = has to be taken into account in order to calculate the 
distribution at higher via DGLAP evolution equations [3]. We show that the origin of 
such low Qo initial distribution can be naturally related to the phenomenon of dynamical 
mass generation of gluons. 

The paper is organized as follows: in the next section we introduce the DGM model and 
show how the semihard gluon-gluon scattering can be properly eikonalized and linked to pp 
and pp diffractive scattering. Our results for h{Q'^) at moderate region are presented 
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in the Sec. Ill, where we also address the question of extrapolating the gluon distribution 
function beyond the infrared range. In Sec. IV we draw our conclusions. 



II. GLUON DISTRIBUTION AND RISING CROSS SECTIONS 

In order to set a natural moderate scale and produce a reliable determination of the 
factor hlQ"^), we explore the small x behavior of g{x^Q'^) = h{Q'^)x~^ in one eikonalized 
QCD- inspired model improved with a dynamical gluon mass |6|. This model (previously 
referred to as DGM model) provides a consistent calculation of total cross sections consonant 
with analyticity and unitarity constraints, and has already been applied successfully to 
describe pp and pp diffractive scattering data 6| , as well as data on 7J9 photoproduction and 



17l |. In the eikonal representation the total cross section is 



hadronic 77 total cross sections 
given by 

atot{s) = An f bdb[l-e-^'^'''^cosxAb,s)l (1) 

where s is the square of the total center-of-mass energy, b is the impact parameter, and 
x{b, s) = Xnib, s) + iXiib, s) is a complex eikonal function. In the DGM model we write the 
even eikonal as the sum of gluon-gluon, quark-gluon, and quark-quark contributions: 

X^{b,s) = Xqq{b, s) + Xqgib, s) + Xgg{b, s) 

= i[agq{s)W{b; fiqq) + aqg{s)W{b; figg) 

+ aggis)Wib;fIgg)], (2) 

where Xpp{.b,s) = x~^{b,s) ± x~{b,s). Here W{b]fi) is the overlap function in the impact 
parameter space and crij{s) are the elementary subprocess cross sections of colliding quarks 
and gluons = q,g). The overlap function is associated with the Fourier transform of a 
dipole form factor, 

W^(6;/x) = ^(Mf i^3(/i&), (3) 
yovr 

where K^i^x) is the modified Bessel function of second kind. The odd eikonal x~{b, s), that 
accounts for the difference between pp and pp channels, is parametrized as 

X-(6,5) = C-S^e-/^iy(6;/z-), (4) 



where rUg is the dynamical gluon mass (which will be discussed afterward) and the param- 
eters C~ and fi~ are constants to be fitted. The factor S is defined as 



(5) 

"a 

with the dynamical coupling constant set at its frozen infrared value. The eikonal func- 
tions Xqqi^: s) and Xgg(^) s), needed to describe the low-energy forward data, are parametrized 
with terms dictated by the Regge phenomenology. The formal expressions of these quanti- 
ties as well as details concerning the analyticity properties of the model amplitudes can be 
seen in Ref. js]. 

In the DGM model the main contribution to the asymptotic behavior of hadron-hadron 
total cross sections comes from semihard gluon-gluon collisions. The gluon eikonal term 
Xgg{b,s) is written as Xgg(^s) = (^ggi^Wib; figg), where 

agg{s)=f drFgg{T)a^^^{s). (6) 

Here Fgg{T) is the convoluted structure function for pair gluon-gluon (r = 0:1X2), and a^^^{s) 
is the complete cross section for the subprocess gg gg, calculated through the dynamical 
perturbation theory (DPT) [20] , where the effective gluon propagator and coupling constant 
enters into the calculation. These effective quantities were obtained as solutions of the 
Schwinger-Dyson equations (SDE) for the gluon propagator and triple gluon vertex via the 
pinch technique Uj, and have been used in many phenomenological calculations that are 
sensible to their infrared finite behavior jl^ . [isl . In this approach the nonperturbative 
running coupling and the functional expression of the dynamical gluon mass M^(g^) are 
given by 

"'^^'^ ^ /?oln [(g2 + 4M2(g2))/A2]' 



1^ 



a 



In 



In'^ 



2N -. -12/11 

(8) 



respectively, where A(= Aqcd) is the QCD scale parameter, q is the four-momentum transfer 
between the incoming and outgoing parton, and Pq = 11 — ^Uf {uf is the number of flavors). 
The dynamical gluon mass rrig has to be found phenomenologically, and a typical value (for 
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A = 300 MeV) is rUg ^ 500 ± 200 MeV [ll|, Q, ll5|, ll6|, ll7| . The expression of M^{q'^) comes 
out as a solution of the gluon propagator SDE and nig is interpreted as the gluon dynamical 
mass, being related to the value of the propagator at = 0. This mass has a dynamical 
origin, which, in principle, should be calculated in terms of the scale A, and in this work it 
will be used as an input parameter. As discussed in detail in the Ref. jQ], the total cross 
section a^^^(s) = I'j"'"''' {da / di) di for the subprocess go — > go is obtained by integrating 
over 4m^ — s < i < 0. In setting these kinematical limits we have neglected the momentum 
behavior in Eqs. ([7j) and ([8]), since the calculation of the hadronic cross section does not 
depend strongly on the specific form of Mg(g^), but more on its infrared value (i.e. the value 
of nig). A straightforward calculation yields 



12s^ — 55m^s^ + 12m^s^ + QQm^s — 8ml ( s — 3m':, 

9 9 9 — 3 In 



(9) 



9 L" ■■"gl \ ■■"g 

In the calculation of the structure function Fggij) = [g<Sig]{T) we adopt a QCD-motivated 
distribution functional form: 

g{x,Q') = h{Q')^-^^P^, (10) 
where the quantity J = 1 + e > 1 controls the asymptotic behavior of a tnt( s). The factor 



(1 — x)^ is introduced in accord with the spectator counting rules at high x 2l|]. In the Regge 
language the quantity J is the intercept of the Pomeron. As we go up in energy (s ^ A^ and 
transferred momenta ~ m^) the proton will be filled up with more and more semihard 
gluons, and we will be probing the small x region. Therefore we see that the behavior of the 
total hadronic cross section at higher energies is driven by only four parameters: m^, J, h 
and figg (see Eqs. (II]), (E]), ([3]) and (ITOll ). The first two are input parameters whereas h and 
figg are fitting ones. As we will see in detail in the next section, h and figg are determined 
carrying out a global fit to high-energy pp and pp scattering data, namely total cross section 
atot and p parameter data. 

In the DGM model the connexion between the dynamical gluon mass and the scale 
Q"^ = — arises as follows: the QCD improved parton model is not applicable only in the 
region in which the transverse momenta Pt of the colliding hadrons is of the same order 
-mass energy, but also in the region in which ^/s ^ pt provided that 
Hence it is possible to investigate the interface between high pt region 



as their center-o: 



Pt > Aqcd 

and Regge limit one by means of perturbative methods. In the DGM model we introduce the 



physical energy threshold s = Aw? for the final state gluons (see expression (ED ), a ssuming 



that these are screened gluons, in a procedure similar to the calculation of Ref. [22|. In this 
way, adopting the relation Px = Q"^, we establish the expected particle production threshold 
Q = 2mg. Note that the gluon distribution function ( |T0|) reduces to g{x,Q'^) = in the 
limit X ^ X , cis expected by dimensional counting rules, and, more importantly, reproduces 
the small x QCD prediction 

gix,Q') = hiQ')x-'. (11) 

In this context, given the high-energy nature of the data to be fitted in the following analysis, 
we are specially concerned with the determination of g{x,Q) in its small x region and 
transferred momenta of the same order as the dynamical gluon masses. As we shall discuss 
in the next section the cross section described by Eq. ([6]) is totally dominated by the small 



X region and by processes at ^ 4m^. However nig ^ 2Aqcd , as discussed in Ref. |ll | 
(although the uncertainty in this value is of 0{Aqcd)), and we are in a moderate region of 
momenta, and above the nonperturbative region determined by the QCD scale (Aqcd)- 



III. THE GLUON DISTRIBUTION AT MODERATE RESOLUTION SCALES 

In the determination of the small scale dependence of the gluon distribution we explore the 
small X behavior of g{x, Q) in one eikonalized QCD-based model improved with a dynamical 
gluon mass. We carry out global fits to high-energy pp and pp scattering data above y/s = 
10 GeV in order to determine phenomenological values for the factor h as a function of 
dynamical gluon mass nig. Note that at this value the dynamics responsible for the 
cross section increase with energy is already effective. These data sets include the total 
cross section {atot) and the ratio of the real to imaginary part of the forward scattering 
amplitude (p parameter). We use the data sets compiled and analyzed by the Particle Data 
Group [2I, with the statistic and systematic errors added in quadrature. In all the fits we use 
a fitting procedure, adopting an interval ~ Xmin corresponding, in the case of normal 
errors, to the projection of the hypersurface containing 90% of probability. In the case 
of 8 fitting parameters (DGM model) this corresponds to the interval ~ Xmin — 13.36. 

As discussed in the previous section, the behavior of the total pp and pp cross sections at 
higher energies is driven specially by four parameters: rrig, J, h and figg. The input values 
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of the nig have been chosen to he in the interval [350, 750] MeV, as suggested by the value 



mn 



mtlllMeV obtained in a previous analysis of the pp and pp channels through the 



DGM model 



61]. This input dynamical gluon mass range is also supported by recent studies 



via the model on survival probability of rapidity gaps in diffractive scattering [16| as well as 

I 

on the photoproduction and the hadronic 77 total cross sections [17]. We stress that if 
we had used rrig as another parameter to be fitted we would reproduce the results of Fig.(l) 
of Ref. j^, indicating that the best gluon mass values would be in the range [300, 750] MeV. 

A consistent input value of J can be obtained from fitting scattering data via Regge 
phenomenology: in the limit of large enough s the gluon-gluon cross section (Jgg{s) obtained 
using the gluon distribution (ITTi) behaves as a Pomeron power law s"'^^: 



lim 



drFgg{T)a^^\s) 



Am'^/s 



(12) 



i.e. the quantity e controls the asymptotic behavior of crgg{s). We see that the asymp- 
totic behavior of <ygg{s) in the DGM model is similar to the one expected from the Regge 
phenomenology, where the quantity J is related to the intercept of the Pomeron, with 
J = 1 + e > 1 From fitting scattering data through an extended Regge model, the 

value of the Pomeron intercept J imposed by the accelerator data currently available is 
J = 1.088±0.011 2^. This value is compatible with the recent result J = 1.085±0.006, ob- 

n 

tained through the use of non-linear trajectories for the meson resonances [25]. Also, a Regge 
fit of the proton structure function data indicates the presence of a second Pomeron, with 
the hard intercept 1 + e = 1.435 [26]. This last value of e is compatible with the one expected 



for the exponent A from the perturbative BFKL approach [27|, namely A = —Os ln2 ^ 0.5. 
However, there is a lack of theoretical precision about the precise value of A once their value 
is dependent upon an infrared cut-off /cq at transverse momenta of gluons 28|, l29| . Moreover, 
it i= abo possible that the phenomenon of dynamical gluon ma. generation p„=h the pertut- 
bative A value to a smaller one as was shown years ago by Ross [30]. More importantly, in 
what follows we will show that the use of a hard value to the Pomeron intercept gives a bad 
fit to atot and p parameter data. It shows that in fact the increase of hadron-hadron total 
cross sections is governed by soft values of J, as already expected from Regge and analytical 



models 



31 



33|. 



The input values for J have been extracted from the Ref. 2J], where an extended 
Regge parametrization was used in order to establish extrema bounds for the soft Pomeron 
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TABLE I: Values of the input parameters and fitting ones of the DGM model resulting from the 



global fit to the pp and pp scattering data. 



m„ fMeVl 


Li„„ [GeVl 




h 


A, / 




J = 


1 + e = 


1.085 




300 


0.741±0.059 




0.162±0.009 


1.112 


350 


0.782±0.064 




0.226±0.013 


1.111 


400 


0.811±0.065 




0.292±0.018 


1.101 


450 


0.857±0.073 




0.357±0.021 


1.114 


500 


0.781±0.061 




0.423±0.028 


1.115 


550 


0.878±0.072 




0.525±0.037 


1.107 


600 


0.920±0.074 




0.634±0.041 


1.108 


650 


0.909±0.071 




0.708±0.052 


1.107 


700 


0.868±0.078 




0.841±0.054 


1.124 


750 


0.956±0.078 




0.993±0.065 


1.124 


400 


J = 

0.679±0.139 


l + e = 


1.435 

0.017±0.001 


1.950 



intercept. In that analysis the intercept of the soft Pomeron were determined exploring 
the systematic uncertainty coming from the discre 3ancy concerning the results for a^f,^^ at 
a/s = 1.8 TeV reported by the CDF Collaboration [3J| and those reported by the E710 ^] 
and the E811 361] Collaborations; the values for the lower and upper bounds are J = 1.085 
and J = 1.095, respectively, and represent the bounds for the soft Pomeron intercept imposed 
by the accelerator and cosmic ray data currently available [24] . We shall restrict ourselves to 
the above range of J values. Values outside this window provide bad fits to the experimental 
data. 

The results of our global fits to the data set discussed above are depicted in Fig. 1. We 
note a fast increase in the h values with the dynamical gluon mass. This behavior does not 
depends on the specific value of the soft intercept, since the results are totally compatible for 
the two intercept bounds, namely J = 1.085 and J = 1.095. The /DOF values obtained 
in the global fits with J = 1.085 are relatively low, as show in Table I (the x^/DOF values 
for the input J = 1.095 are similar). These results (for 148 degrees of freedom) indicate 



9 



the excellence of the fits and show that it is quite unlikely that the soft Pomeron does 
not dominate the behavior of total cross sections at least in the energy region that we 
are considering. In other words, our fits do not support a hard Pomeron, as indicated in 
the Table I: the use of the hard input value J = 1.435 results in a significantly worse fit 
{x^/DOF = 1.95) and shows that in this case the DGM model does not accommodate the 
data set used in the fitting procedure. Moreover, the input choice J = 1.435 results in a 
very low value for the factor h, namely h = 0.017 ±0.001. As we will see in the next section, 
this value is approximately 1/10 of the expected one [37|. 

The small x behavior of xg{x,Q^) = h{Q'^)x~'^ at = 1 GeV^, taking into account our 
h result for Trig = 500 MeV as indicated in Table I, is given by 

xg{x) = (0.423 ± 0.028)x-(°°^^±°-°^), (13) 

where we have adopted the intercept bound 1 + e = 1.085^ which is consistent with the recent 
value of J obtained from spectroscopy data analysis [25]. This behavior can be compared 
with the MRST result for the LO parametrization of g{x) at = 1 GeV^ at very small x 



xg^'^'^'^ix) ^ 3.08x°-i°; (14) 

at x' = 10^^ (the kinematical limit of the MRST distributions) our distribution gives 
x'g{x') = 1.125]'^o'}33, whereas x'g'^^^^^{x') = 0.974. This numerical agreement reinforces 
the assertion that the hadronic scattering is driven mainly by soft Pomerons at scales ^ 1 
GeV^. This is not the case for higher scales, as we can see in Figure 2, where we compare our 



gluon distribution with the MRST 



39| leading-order gluon 



and CTEQ6L1, CTEQ6L 
distributions at = 1.69 GeV^ and = A GeV^. These parton distributions are evolved 
with LO splitting functions using DGLAP equations from the initial scales Qq = 1 GeV^ 
(MRST) and Ql = 1.69 GeV^ (CTEQ6), and are obtained using tree-level formulas for the 
hard cross sections. The CTEQ6 (MRST) parton sets are valid in the kinematical intervals 
10-6 < X < 1 (10-5 < a; < 1) and 1.3 GeV < Q < 10^ GeV (1 GeV < Q < 10^ GeV). 
The CTEQ6L distribution uses the same asiQ"^) coupling as the standard next-to- leading 
order (NLO) CTEQ fits with a,(M|)=0.118 whereas CTEQ6L1 one uses the LO formula for 
a,(Q2) with Aq^cd"^ = 0.215 GeV (with corresponds to as(M|) = 0.130). The LO MRST 
distributions uses a value of AjS^"*^ = 220 MeV corresponding to a.(M|) ^ 0.130. In the 
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DGM model we adopt the QCD parameter value A(n/ = 4) = 300 MeV, which corresponds 
to as{Ml) ^ 0.130. The results depicted in Fi gure 2 indicate that the small-x gluons at 
higher scales are mainly generated radiatively: at the value = 4 GeV^ the DGM 
gluon distribution and the CTEQ/MRST ones are discrepant at greater than a factor of 5 
at X = 10-^ ahhough the DGM and MRST results are compatible at = 1 GeV^ scale. 
The DGM gluon distributions at higher scales can be calculated by fitting the h results 
to the functional for h{Q'^) = A + BQ^, as shown in Figure 3, where the best fit gives 
A = 0.121 X 10"^ and B = 0.170 x 10"^ MeV"^ jg ^ear that the functional form of 
h{Q'^) is not known and we stress that the choice we made above is just a guess, as well as 
it should be kept in mind that the values that we are using in this function are of the 
order of 4m^. 

It is worth mentioning that mass values nig < 500 MeV are favored in other calculations 
of strongly interacting processes, as we will discuss later. In this way, we consider the 
more stringent moderate-Q^ interval 0.36 GeV^ < ^ 1 GeV^ (corresponding to the 
gluon mass interval 300 MeV < rrig < 500 MeV) as a more reliable one for the DGM 
gluon distribution. This new threshold interval for the production of dynamically generated 
gluons is corroborated by our results at > 1.69 GeV^ {nig > 650 MeV) since our gluon 
distribution is smaller than the LO CTEQ6 and MRST ones. It is a clear indication that in 
the DGM model threshold values nig > 500 MeV inhibit the production of nonperturbative 
gluons. 

At even smaller scales, for instance < 0.36 GeV^, the x behavior of the gluon dis- 
tribution is due to the dynamics of the bound state proton and is an unsolved problem of 
nonperturbative QCD. Hence we guess a behavior for g{x, Q^) that provides a picture more 
consistent with the expected behavior of the nonperturbative dynamics of the bound state 
nucleon 40|, |4l|]: we assume the frozen form g{x, rf) for momenta Q < rj = 300 MeV. In this 



way, we can write down a phenomenologically useful expression to the gluon distribution 
function in the range Q < 1 GeV: 

gix, Q') = h{v') 9iv' - Q') + h{Q') e{Q^ - r^'). (15) 

This phenomenological gluon distribution takes into account the fact that at small 
scales {Q < 300 MeV) the usual methods of analysis of parton distribution based on per- 
turbative QCD are not applicable. It is clear that low-Q^ photons do not resolve partons 
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in the proton and the deep inelastic scattering (DIS) formahsm cannot be apphed straight- 
forwardly, and probably a generalized vector dominance model, as discussed by Alwall and 



Ingelman 42(], is more appropriate at low . Generally speaking we could say that we 
should expect a smaller gluon distribution function at moderate than the one obtained 
through a strictly perturbative technique (the CTEQ approach, for instance), because of the 
propagator and coupling constant softening in the infrared which will appear in the cross 



section calculation 
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IV. CONCLUSIONS 

From the diffractive formalism underlying QCD-based models we have showed how the 
rise of total hadronic cross sections can be naturally connected to the small x behavior of the 
QCD-predicted gluon distribution function g{x,Q) at moderate Q'^. We have performed a 
precise determination of the h^Q"^) factor that appears in the QCD-predicted gluon distribu- 
tion in the framework of the DGM eikonal model. This QCD-based model is instrumental for 
a reliable determination of the small x behavior of g{x^ Q) at moderate since diffractive 
processes involve both semihard and soft scales resulting in a complicated interplay between 
perturbative and nonperturbative effects. More specifically, we have determined the behav- 
ior of the gluon distribution g{x,Q'^) = h{Q'^)x~'^ in the moderate Q interval Q ^ 1 GeV. 
The connexion between the resolution scale Q and the gluon dynamical mass rrig has been 
estabhshed by the particle production threshold Q = 2mg. 

Our results show that at momenta ^ 1 GeV^ the small-x gluons are mainly generated 
radiatively. The coefficient h is dependent as discussed at the end of the last section. 
We verified that J values like J = 1 + e ^ 1.08 are preferred instead of the large (hard 
perturbative Pomeron) value J = 1 + e ~ 1.43. However the rapid increase of h{Q'^) with 
the momentum is not enough to produce a number of nonperturbative gluons as high as 
the perturbative ones. For this is probably necessary a rapid increase of the J values with 
the momentum in such a way that soft values are preferred at moderate and hard ones 
would be necessary only at high . We call attention to the fact that this statement is 
corroborated by a MRST group analysis of parton distributions of proton 43| . From fitting 
the sea quark (S) and gluon (G) distributions of the default MRST partons to the forms 
xfi{x, Q^) = A{Q)x~^^^'^^^ as X ^ 0, i = S,G, they have look at how the exponents and 
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Xg vary with Q^. It was observed that as increases from the input scale Qo = 1 GeV^ 
the valence-hke character of the gluon rapidly disappears due to evolution being driven by 
the much steeper sea. For higher values of the gluon exponent Xq increases rapidly and 
becomes higher than the sea quark exponent A5, since the gluon drives the sea quark via the 
g ^ qq transition. By taking only positive values for Ac we observe that it starts at a value 
where the gluon distribution is almost flat, that is Ac ~ 0, and by ^ 4 GeV^ it has the 
value Xg = 0.2. The hard value Xg = 0.435 is achieved at Q'^ ^ 700 GeV^. Therefore our 
results are consistent with the conventional view that the steeply-rising gluon component is 
absent at moderate and as increases it is generated through pQCD evolution. 

The small x behavior of the gluon distribution (specifically the value of h) was performed 
in a different approach by Cudell, Donnachie and Landshoff (CDL) [37], within the Land- 
shoff and Nachtmann Pomeron model (LN) {44, 
distribution function is given by 



461 ]. In this model the small x gluon 



xg(x) = h = ^ [ dk^k^D^-e), (16) 
Stt J 

where the gluon propagator (D) is the effective one with the perturbative part subtracted 
from the integral. Notice that Eq. f|T6l) does not contain the x"*^ contribution which has to 
be put by hand in the LN model. This calculation with Cornwall's propagator, where we 



can substitute a 
CDL approach 



)y as, gives h 0.15: this value is basically equal to the value obtained in 
371]. A comparison of our results with the CDL one shows that our values of 
h are quite reasonable. More explicitly, shows that the CDL result of h is compatible with 
the DGM one at mass scales on the order of rn^ ~ 300 MeV (see Table I). Dynamical gluon 
masses of that order where obtained in Ref. [161] in order to obtain the survival probability 
(|S'p) of rapidity gaps in diffractive scattering from production of Higgs boson via WW H 
fusion processes. These survival factors are compatible with recent results of (jS'p) obtained 
by Khoze, Martin and Ryskin through a two-channel eikonal model which embodies pion- 
oop insertions in the Pomeron trajectory, diffractive dissociation, and rescattering effects 



471 ]. In addition, in the analysis via DGM model of the 7J9 photoproduction and the hadronic 
77 total cross sections, both derived from the pp and pp scattering amplitudes assuming 
vector meson dominance, there is room for smaller values of rrig beside the adopted mass 



scale TJig = 400 GeV IT]]. All these results suggest a most restrictive mass range than the 



one we have considered up to now. Given the behavior of g{x, Q^) obtained through the 

13 



DGM model, we can be confident in the mass range mg < 500 MeV as a more reliable one. 
This corresponds to the scale interval "^1 GeV^. It is clear that the dynamical gluon 
mass is a concept that has not been made totally compatible with perturbative QCD, this 
will probably demand time and a great effort, although the introduction of this concept is 
clearly consistent with the experimental data and indicate a more appealing view to the 
behavior of g{x) at small x and moderate Q^. 

In summary, we think that the DGM model is an useful tool in determining the behavior 
of the gluon distribution function at moderate by merging in a natural way the soft 
and the semihard components of hadronic scattering processes. We emphasize that it is 
a nontrivial result that our gluon distribution behavior at small x is on the order of the 
MRST one at = 1 GeV^. This is achieved with only three parameters, namely m^, J 
and h. The natural relation between the momentum scale Q and the dynamical gluon mass 
rUg argues the dynamical mass generation of partons to be the main mechanism behind the 
nonperturbative dynamics of the valence-like parton distributions present at low resolution 
scales. 
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FIG. 1: The factor h as a function of the gluon dynamical mass rUg determined by means of the 
DGM model. 
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FIG. 2: The gluon distributions at = 1.69 and 4.0 GeV^ from the DGM model (sohd curves). 
The dotted, dashed and dotted-dashed curves are the distributions corresponding to the MRST, 
CTEQ6L1 and CTEQ6L gluon sets, respectively. 
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FIG. 3: The factor h ss a function of the gluon dynamical mass irig. The sohd curve is obtained 
through the function h{Q^) =A + BQ^, with A = 0.121 x 10"^ and B = 0.170 x 10"^ MeV-^. 
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